Collision number statistics for transport processes 
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Many physical observables can be represented as a particle spending some random time within a 
given domain. For a broad class of transport-dominated processes, we detail how it is possible to 
express the moments of the number of particle collisions in an arbitrary volume in terms of repeated 
convolutions of the ensemble equilibrium distribution. This approach is shown to generalize the 
celebrated Kac formula for the moments of residence times, which is recovered in the diffusion limit. 
Some practical applications are illustrated for bounded, unbounded and absorbing domains. 



The evolution of complex physical systems can of- 
ten be described in terms of particles undergoing ran- 
dom displacements. Randomness may reflect the in- 
trinsic stochastic nature of the transport process, or re- 
sult from the uncertainty affecting the interactions be- 
tween the travelling particle and the surrounding envi- 
ronment [H, Q. In this context, assessing the statistics of 
the residence time iy spent by the random walker inside a 
given domain V plays a key role in many practical prob- 
lems, encompassing areas as diverse as research strate- 
gies, market evolution, percolation through porous me- 
dia, and DNA translocation through nanopores, to name 
only a few This has motivated a large number 

of theoretical investigations over the last decade, cover- 
ing both homogeneous and heterogeneous, scale- invariant 
media 0-0] ■ In the former case, the dynamics of the 
walker is typically modelled by regular Brownian mo- 
tion [lfj, whereas in the latter one resorts to anomalous 
diffusion [l"0 ] . For Brownian motion, in particular, a sem- 
inal work developed by Kac [IH has allowed all the mo- 
ments of the residence times to be evaluated by resorting 
to convolution integrals over the ensemble equilibrium 
distribution of the walkers, for arbitrary boundary condi- 
tions on V [13]. When the particle is lost upon touching 
the boundary dV of V, the residence time is then usu- 
ally called first-passage time Q, and this quantity has 
been extensively studied for both Brownian motion and 
anomalous diffusion || [^, [13] ■ 

However, in many realistic situations, stochastic trans- 
port is dominated by finite-speed effects, so that the dif- 
fusion limit is possibly not attained, and the dynamics of 
the walker is better described in terms of the Boltzmann 
equation, rather than the (anomalous) Fokker-Planck 
equation [ijj Examples are widespread, and arise 
in, e.g., gas dynamics, neutronics and radiative transfer, 
electronics, and biology fl^ - l20l ]. In all such systems, the 
stochastic path can be thought of as a series of straight- 
line flights, separated by random collisions, as in Fig. [1] 
A natural variable for describing the walker evolution is 
therefore the number of collisions riy within the observed 
volume. Application of the diffusion approximation to 
the characterization of the residence times distribution, 
which amounts to assuming a large number of collisions 
in V, might lead to inaccurate results [2l|. In the present 



Letter, we address the issue of generalizing Kac approach 
to random walkers obeying the Boltzmann equation. We 
derive an explicit formula relating the moments of ny to 
the equilibrium distribution of the walkers, for both scat- 
tering and absorbing media. Knowledge of higher order 
moments allows on one hand estimating the uncertainty 
on the average residence time, and on the other hand re- 
constructing the full distribution of the collision number. 
We show that when V is large as compared to the typical 
size of a flight, so that the diffusion limit is reached, Kac 
formula is recovered. 

Methodology. Consider the random walk of a particle 
starting from a point-source located in ro. At each colli- 
sion, the particle can be either scattered, with probability 
p, or absorbed (in which case the trajectory terminates). 
For the sake of simplicity, we assume that scattering is 
isotropic. We denote by r the position of the walker 
entering a collision, as customary. Let n(r, r') be the 
probability density of performing a displacement from r' 
to r, between any two collisions. Define the transport 
operator 7r[/](r) 



7r[/](r)= f 7r(r,r')/(r')dr', 
Jy 



(1) 



over a d— dimensional volume V. We can then express 
the propagator \t(r, n|ro), i.e., the probability density of 
finding a particle in r at the n— th collision (starting from 
ro), as 



*(r,n|r ) = p 



n-l^nr 



(r,r ), 



(2) 




FIG. 1: A random walk starting from ro and performing a 
limited number n\> of collisions in a region V (with transparent 
boundaries), before being absorbed at r„. 
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FIG. 2: Moments (n v )(ro) (stars) and (riy)(ro) (circles) for 
3d a — 2 gamma nights with scattering probability p = 1 in 
a volume V with transparent boundaries. Theoretical predic- 
tions: solid lines; Monte Carlo simulations: symbols. 
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FIG. 3: Moments (n\>)(ro) (stars) and (riy}(ro) (circles) for 
Id exponential flights with scattering probability p — 0.5 in 
a volume V with transparent boundaries. Theoretical predic- 
tions: solid lines; Monte Carlo simulations: symbols. 



where ir n [/] (r) is the n-th iterated operator 

7T n [/](r)= / [ ir(r,r n )...ir(r 2 ,r 1 )f(r 1 )dr 1 ...dr n , (3) 
Jv Jv 

and S = S(r — ro) is a short-hand for the initial point- 
source condition. 

The probability of performing ny collisions in the vol- 
ume V is related to the propagator by 

V(n v \r Q ) = / dr^(r,n v \r )- dr*(r, n v + l|r ). (4) 
Jv Jv 



The moments 



+oo 



l )(r ) = ^r(n v \r ) 



(5) 



depend on the boundary conditions on dV, which affect 
the functional form of the propagator. The absence of 
boundary conditions corresponds to defining a fictitious 
('transparent') volume V, where particles can indefinitely 
cross dV back and forth. On the contrary, the use of 
leakage boundary conditions leads to the formulation of 
first-passage problems [7|. 

We introduce now the collision density ^(r|r ) 



N 



*(r|r ) = Jim V^(r,n|r ), 



(6) 



which intuitively represents the equilibrium (stationary) 
particle distribution [l6j. From the definition of \P(r|ro), 
it follows immediately that 



(4)(r ) = f dr*(r|r ), 
Jv 



(7) 



i.e., the integral of the collision density over a volume V 
gives the mean number of collisions within that domain, 
hence the name given to \&(r|ro). 

Higher order moments of rty can be obtained as follows. 
Define the operator 



vf[/](r) = / tf(r|r')/(r')dr'. 
Jv 



By making use of the Neumann series 

oo 

£p-V[/](r) = — [/](r), 



n=l 

we have then 



*[/](r) 



1 — piT 



[/](r) 



(8) 



(9) 



(10) 



and in particular ^(r|ro) = *&[5](r, ro). Now, combining 
Eqs. J5J and (gj), we get 



1 

P Jv 



drLi_ m (p7r)(l-p7r)[J](r,r ), (11) 



where Li„(.r) = J2kLi % /k s is the polylogarithm func- 
tion [2?!]. When to is a non-negative integer, the polylog- 
arithm is a rational function, namely, 



Li 



i- m (x) = ^2 k\s rn+1 , k+1 ( - 



fc=0 

where the coefficients 



fc+i 



1 k 



(k-iy 



(12) 



(13) 
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FIG. 4: Moments (ny)(ro) (stars) and (ny)(ro) (circles) for 
Id exponential nights with scattering probability p = 1 in a 
volume V with leakage boundaries. Theoretical predictions: 
solid lines; Monte Carlo simulations: symbols. 



are the Stirling numbers of second kind [22J. Thanks 
to the recurrence properties of the Stirling numbers, 
Eq. (fTTj) gives 

m 

«)(r ) = - / rfrVfe! Sm)fc /^[5](r )r o). (14) 

We introduce then the repeated Kac integrals 

C k (r Q ) = k\ [ dr k ... f dr 1 *(r jfc |r fc _i)...*(ri|r d ), (15) 
Jv Jv 

which arc dchned as fc-fold convolutions of the collision 
density \l/(r|ro) with itself 12, It follows the equiv- 
alence 



C k (r ) = k\ / dr* fc [<J](r,r ), 
Jv 



(16) 



with C fc+ i(r ) = fc*[C fc ](r ) and C (r ) = l ro ev, 1 being 
the characteristic function. The convergence of the in- 
tegrals Cfc(r ) depends on the features of the underlying 
stochastic process as well as on boundary conditions. For 
instance, the persistence property of walks in d < 2 (Toj 
implies diverging Cfc(r ) for transparent V in absence of 
absorption 13]. We finally obtain the central result of 
this Letter, i.e., the explicit formula for the moments of 
the collision number in V 



W)(ro) = -f; Sm ^(r„). 



(17) 



Thanks to linearity, Eq. (fTTj) allows expressing 
(riy)(ro) as a combination of m Kac integrals, k = 
l,...,m, each given from Eq. (|15[) . In particular, for 



rn = 1 we recover Eq. ([7]), since Si,! 
for m = 2, S2.1 = S22 = 1, so that 



1. Furthermore, 



(nl)=p2\ f dr 2 f dri*(r 2 |r 1 )*(r 1 ,r ) + <nt>. (18) 
Jv Jv 

Let G(z\rp) be the moment generating function of 
■p(nv|i"o) [10j . By definition we have the moment ex- 
pansion 



G(z\r )=J2(n^(r ) 



ml 



(19) 



The small-z expansion of G(z|ro) reads then G(z\tq) ~ 
1 + (nh)(To)z, which for the Tauberian thorems corre- 
sponds to the large-ny behavior. It follows the exponen- 
tial tail 



T(n v \r ) ~ e ~ nv/{n ^ ){ra \ 



(20) 



provided that Ci(ro) is finite. 

Diffusion limit and Kac formula. Suppose for the sake 
of simplicity that the walker moves at constant speed v, 
so that the time spent between any two collisions is tj = 
|r.j — Ti_i\/v. For isotropic walks, 7r(r, r') = tt(£ — \r— r'|) 
thanks to the spherical symmetry. Then, flight times 
are identically distributed, and obey U ~ w(U), where 
w(ti) = Cl d J ^ 1 7r(£)5(i i - t/v)d£, n d being the surface 
of the unit sphere. The diffusion limit of the transport 
process described above is obtained by letting the typical 
flight length a (the standard deviation of jump sizes, as 
customary) and the average inter-collision time r = (ij) 
shrink to zero in such a way that the ratio D — o~ 2 /t 
converges to a constant, namely, the diffusion coefficient. 
When t and a vanish, the collision number in V diverges, 
whereas the quantity 



"V 

i=l 



(21) 



converges to the residence time in the volume. Actu- 
ally, iy should take into account also additional terms 
due to boundary conditions. However, as r — ¥ and 
iiy — > 00, the trajectory will almost surely have a turn- 
ing point touching the boundary, so that corrections to 
Eq. (|2"T1) can be safely neglected. Let now Q(iv| r o) 
be the distribution of the residence times. Under the 
previous assumptions, in the Laplace space we have 
Q(s\r ) = J exp(-st v )Q(t v \r )dt v = w(s) n v . Any arbi- 
trary w(s) with finite r has an expansion w(s) ~ 1 — sr 
when t — > 0. Then we have Q(s|rrj) — e~ nyST , which 
implies Q(£y| r o) ~ S(t\> — nyr) for small r. It follows 
that 



<*v>M 



■ l )(ro), 



(22) 



with (n™ )(r ) given by Eq. (fT7|> . If we rescale the space 
variable r by cr, each of the terms of the sum in Eq. (|17[) 



4 



carries a contribution a 



-2k 



In the diffusion limit 



therefore recover the celebrated Kac formula [12|, [1 
the moments of the residence times of Brownian motion 



we 
for 



<C)(ro) 



Cm(rp) 
D m 



(23) 



because all other terms in the sum vanish when t — > 
and a — > 0, and s m>m = 1. Moreover, we have the re- 
cursion property D{t™ +1 )(r ) = m*[(t^)](r ), starting 
from (ty)(r ) = l ro ev- As for the residence time distri- 
bution, from Eq. ([2TJ1) at long times we have Q(iv| r o) — 
exp(-I>tv/Ci(ro)). 

Discussion and perspectives. For simple geometries 
and W(r|ro), formula (|17p is amenable to analytical solu- 
tions. Here we illustrate some significant cases, V being 
a d-sphere of radius R centered in 0. Isotropic gamma 
flights with kernel n(£) = £ a - d exp(-£)/n d T(a), a > 0, 
are a widespread transport process and describe, among 
others, search strategies in biology (23[. When a = 2, the 
3d scattering collision density assumes a simple form, 



(r|r ) = 



1 



47r|r - r 



It follows 



(4>M 




and 



<«v>(ro) 



' 25fl 4 -10-R 2 rg+rg 



r <R 
r >R 



+ {nl)(r ) r <R 



15 r 



r >R 



(24) 



(25) 



(26) 



Comparisons with Monte Carlo simulations with 10 5 par- 
ticles are shown in Fig. [2] Exponential flights (a = 
1) arise when the scattering centers are uniform, so 
that inter-collision distances obey a Poisson distribution. 
Such process is crucial for understanding, e.g., radia- 
tion propagation [lfjj]. The transport kernel reads 
ir(£) = £ d exp(— £)/fld- In Id systems, the collision 
density for absorbing V and transparent boundaries reads 



*(r|r ) 



271— p ' 

whereas for scattering V and leakage boundaries 



*(r|r ) 



(27) 



(28) 



R+l 24]. The moments (rc^)(r ) are again [22 



with r L 

easily obtained from Eq. (fTTjl . The resulting formulas 
are rather cumbersome and will not be presented here. 
Comparisons of exact results and Monte Carlo simula- 
tions with 10 5 particles are shown in Figs. [3] and [4] 



The integrals appearing in Eq. ([i~7|) can be carried out 
numerically for arbitrary complex volumes and \P(r|ro). 
Moreover, the isotropy hypothesis can be possibly relaxed 
by replacing r with a state variable y = {r, 6} account- 
ing for the scattering angle. In this respect, the approach 
presented in this Letter is general, and applies to a broad 
class of transport processes and boundary conditions. Fi- 
nally, we observe that Eq. (fl~7| could be used as an esti- 
mator to infer the properties of the underlying stochastic 
path, which is often not accessible, on the basis of the 
collision statistics in a measure volume. 
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